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ABSTRACT 



In this paper, we show that if the broad line region clouds are in approximate energy 
equipartition between the magnetic field and gravity, as hypothesized by Rees, there 
will be a significant effect on the shape and smoothness of broad emission line profiles 
in active galactic nuclei. Line widths of contributing clouds or flow elements are much 
wider than their thermal widths, due to the presence of non-dissipative MHD waves, 
and their collective contribution produce emission line profiles broader and smoother 
than would be expected if a magnetic field were not present. As an illustration, a 
simple model of isotropically emitting clouds, normally distributed in velocity, is used 
to show that smoothness can be achieved for less than ~ 8 x 10 4 clouds and may even 
be as low as a few hundred. We conclude that magnetic confinement has far reaching 
consequences for observing and modeling active galactic nuclei. 
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1 INTRODUCTION 

Quasar emission lines are important probes of the physics 
of active galactic nuclei (AGN). The primary assumption 
made in spectral line synthesis studies of AGN is that the 
width of a component contributing to the total line profile is 
thermal. The assumption is made regardless of whether the 
component is a single cloud in a discrete ensemble of clouds 
or a differential volume element in a continuous flow. 

Line emission originates in matter at ~ 10 4 K and cor- 
responds to a thermal width of ~ 10 km s _1 . Since broad 
line region (BLR) emission lines have much larger widths 
(FWHM ~ 10 3 to ~ 10 4 km s" 1 ) the width of an individ- 
ual component is often assumed to be negligible compared 
to the total line profile. Based on this assumption and the 
observed smoothness of broad emission line profiles, cloud 
numbers in excess of 10 7 to 10 8 are inferred (Arav et al. 
1997, Arav et al. 1998, Dietrich et al. 1999). 

This need not be the case however. If a magnetic field 
is present, the line width of a contributing element will be 
broadened requiring fewer clouds to produce the observed 
profile smoothness. 

In nature a magnetic field is usually associated with 
non-dissipative MHD waves in energy equipartition with the 
magnetic field. Thus 
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where B 2 /8tt and l/2p<r B are the magnetic pressure and 



MHD wave energy density respectively, p is the mass density 
and ob is the resulting velocity width of the gas. (Arons 
and Max 1974, Meyers and Goodman 1988a, Meyers and 
Goodman 1988b). 

Rees (1987) suggested that BLR clouds are magneti- 
cally confined. Assuming 
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and solving for B gives 
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where nio is the density in units of 10 10 cm" 3 and T4 is the 
temperature in units of 10 4 K. From this we see that only 
a few Gauss are required for confinement. Substitution of 
Equation 3 into Equation 1 and solving for ctb gives a lower 
bound for the line width of magnetically confined BLR gas. 
Thus 
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where tua is one atomic mass unit and Z is the mean atomic 
weight of the gas which, assuming cosmic abundances, is 
taken to be Z « 1.4. This is comparable to the thermal 
width of hydrogen, but is roughly 3.5 times larger than 
the thermal width of carbon (~ 3.2 km s _1 ), and 7.5 times 
larger than the thermal width of iron (~ 1.5 km s _1 ) at 
T4 « 1.0. Thus, even for the minimal confining magnetic 
field there will be significant effects on line transfer. 
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The broadening is considerably amplified however if the 
magnetic field is in equipartition with the gravitational en- 
ergy density so that 

n 

8?r 



GMp 



R 



(5) 



(Blandford and Payne 1982, Rees 1987, Emmering et al. 
1992, Konigl and Kartje 1994, Bottorff ct al. 1997). Al- 
though there is no fundamental reason why Equation 5 
should hold, this equipartition does occur in many environ- 
ments (Rees (1987), Meyers and Goodman (1988a,b)). Here 
R is the radial distance of gas, with mass density p, from 
a central black hole of mass M. In this case, the line width 
will be greater than the virial width, the local velocity field 
will be highly supersonic and MHD broadening may actually 
account for the full width of an emission line. 



2 A SIMPLE MODEL 

To illustrate the effect of non-dissipative MHD wave line 
broadening and smoothing on an emission line profile, we 
consider the extreme case of BLR emission arising from a 
discrete set of identical clouds. This example has immediate 
applicability in the search for discrete clouds or extended 
cloud structures in BLR line profiles. (Arav et al. 1997, Arav 
ct al. 1998, Dietrich et al. 1999) 



2.1 A Discrete Distribution of Emitters 

A simple one dimensional outflow model is constructed in 
which clouds move only along the line of sight. We choose a 
cloud bulk velocity field given by 
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Requiring consistency with Equation 1 and Equation 5 gives 
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Thus the dispersion of an emitting element equals the sys- 
temic velocity at any given radius. (Note: A magnetic ex- 
ample in which va ~ y/2GM/R is the MHD wind model of 
Blandford and Payne (1982). In that paper va is also given 
by vg ~ y B 2 /4:TYp. An association of <jb with vg, however, 
is not pursued.) 

In the model being considered here identical clouds are 
placed randomly in velocity space according to a gaussian 
distribution given by 

1 / vg \ 2 

f (vg Iog ) oc exp-- — (8) 
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where <jg is the dispersion of the cloud distribution in veloc- 
ity space. We loosely associate era with the mass M and an 
emission weighted radius Rq (e.g. the reverberation radius) 

giving o g ~ y // 2GM/RG- Thus og does not include the ef- 
fect of magnetic broadening which must be added separately 
to each cloud. 

We will predict the profile of the 1549 A line of C IV. 
The surface emissivity of a cloud, e(C iv), is given by an 
analytical fit of C IV emission for an amalgam of BLR cloud 
densities as prescribed in Baldwin et al. (1995). The fit is 




Figure 1. Plot of the distribution /(vg/ug) normalized to 
/(0.0) = 1.0 and including a (vGfo'G)cut = 0.24 cutoff ((a), 
thick line). Also shown is g(— 0.32, vq/ctg) and g(2.05, vg/<?g) 
((b) and (c) respectively, thin lines). Both are normalized so that 
g(2.05,0.0) = 1.0 



given in Bottorff et al. (1997) and is reproduced here for 
convenience. 
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Here $ 18 (H) = $(//)/10 18 where $(#) is the hydrogen 
ionizing photon number flux in cm -2 s _1 . Following Netzer 
and Laor (1993) we assume that lines are suppressed with 
the onset of grain formation at &i$(H) — 1.0 so we assign 
zero emissivity to clouds exposed to this flux or less. This is 
satisfied for 
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where Rg,io is the radius Rg, written in units of 10 
light days, L45 is the bolometric luminosity in units of 
10 4B erg s _1 , and (vg/og)^ is the value of \vg/ug\ be- 
low which the emissivity is defined to be zero. We take 
Rg,io ~ 10 and L45 « 0.27 (values corresponding to 
the' Seyfert 1 galaxy NGC 5548, Bottorff et al. 1997) so 
(vg /era) cut ~ 0.24. For comparison, an example of a quasar 
is 3C390.3 which has R G ,w ~ 6.3 and L45 ~ 1.8 (Wamsteker 
et al. 1997) giving (vg /<?G)cut ~ 0.43. The cutoff in emis- 
sivity is equivalent to truncating the distribution f(vG I &g)- 
The thick curve in Figure 1 shows /(vg/ctg) normalized to 
/(0) = 1.0 and truncated for \v G /o G \ < 0.24. 

Our simulation used a total of 84,000 clouds, estimated 
from typical BLR cloud column densities, particle densities, 
and the size and covering fraction of the BLR. To simulate 
MHD wave broadening, each cloud was given a gaussian line 
profile centered at a randomly selected value of vg/og de- 
noted as Vi/oG and assigned a dispersion equal in magnitude 
to Vi /og so as to be consistent with Equation 7. The cloud 
line profile, g(vi/aG,VG /o~g), is thus 
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A cloud with small \vi/<JG\ (but still larger than 
(vg /cG^cut), has a relatively narrow width and makes 
a smaller small flux contribution to the total line pro- 
file as compared to clouds with larger |v;/cxg|. For com- 
parison Figure 1 also shows two cloud profiles, namely 
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Figure 2. Simulated BLR portion of an emission line profile for 
our simple 1-D model. For this simulation 84,000 clouds were used 
with (v G /o- G ) C ut = 0.24. 

g(— 0.32, vg/og) and g(2. 05, vg/o'g)- Both have been nor- 
malized to g(2.05, 0.0) = 1.0. It is apparent in the figure 
that, individually, clouds with high |vi/(Tg| outshine those 
with low \vi/a G \- For the two cloud profiles shown in Fig- 
ure 1 the larger \vi/a G \ cloud is 10 times more luminous. 
On the other hand there are many more low velocity clouds 
than high velocity clouds due to the distribution /(vg/og)- 
The resulting profile, due to the accumulation of all 84,000 
clouds is shown in Figure 2 though only 81 percent of the 
clouds actually contribute the line profile due to the cut- 
off. The Profile represents the BLR contribution of an AGN 
emission line. 

Analysis of this model shows that lowering the number 
of clouds to 300 has little effect on the line profile though it 
does become somewhat less symmetric due to the sensitiv- 
ity of the profile to large individual contributions from high 
velocity clouds. The effect of MHD wave broadening is ap- 
parent when the full width at half maximum (FWHM ) of the 
profile in Figure 2 is compared to the FWHM of /(vg/og)- 
The line profile is 1.6 times wider than f (vg/o'g)- The same 
model with 84,000 clouds but using (v G /cr G )cut = 0.43 (e.g., 
parameters for 3C390.3) yields a similar profile to the one 
shown in Figure 2. It is somewhat wider (the FWHM is 1.7 
times larger than the FWHM of f (vg/og)) but has a flat 
plateau nearly twice as wide. 

In the following section we interpret the above results 
in terms of current AGN research and make suggestions for 
future avenues of study. 

2.2 Can the Cloud Nature be Determined? 

This simple model brings a whole series of issues to the 
study of AGN phenomena. Attempts to infer cloud num- 
bers from broad line profiles need to include the possibility 
of MHD wave line broadening before making conclusions 
about whether BLR clouds are continuous or discrete. If in- 
dividual clouds are sought, the search needs to be redirected 
away from the line wings and toward the line core since, 
by Equation 4, the line broadening, <tb, will be weaker at 
smaller values of v. 

An alternative approach to AGN cloud counting could 
be to use a principle component analysis fitting approach. 
The gaussian basis functions used to fit a profile can be given 



widths proportional to their offset from zero systemic veloc- 
ity. The minimum number of gaussians required to find an 
acceptable fit to the line profiles will be an estimate of the 
minimum cloud number. Components may be tested in AGN 
that have had extensive spectral monitoring, e.g. NGC 5548, 
NGC 4151, 3C390.3 and 3C273. A sequence of spectra, cov- 
ering a time span shorter than the BLR dynamical crossing 
time but longer than a characteristic continuum variability 
timescale, can be fit. If the component number and locations 
in the sequence do not significantly vary, then that would 
be evidence in favor of the components being actual clouds. 
The ensemble could be further tested against various kine- 
matic models by tracking detected clouds over a few crossing 
times. There is already evidence to suggest that line profiles 
of clouds are a complex amalgam of emitters. Multiple com- 
ponents are often required to fit line profiles. For example, 
four or five components are required to fit the H/3 line of 
Ark 120 (Korista 1992) and H/3 in NGC 5548 shows three 
seemingly independent time variable components (Wanders 
and Peterson 1996). We note, for clarity however, that it 
is not suggested that the components presented in those 
papers necessarily represent actual clouds since the fitting 
algorithms are designed for efficiency and have no physical 
meaning. We do however, wish to emphasize the potential 
for future reanalysis of available data. 

With regard to magnetic broadening and long term pro- 
file variability, consider the "shoulders" of the H/3 line of 
NGC 5548 reported by Wanders and Peterson (1996). They 
note that shoulders "do not appear to move systematically 
in radial velocity but appear to come and go at approx- 
imately fixed wavelengths." In terms of our simple model 
this behavior can be explained by the movement of a few 
relatively high velocity clouds. Since both the emission and 
line width is large in these clouds we would expect that the 
wings vary on a time scale of the order of the BLR crossing 
time, which they estimate to be ~ 1.6 yr. This is indeed the 
case. The line core, in our model, however, is dominated by 
many dimmer clouds. Thus the profile core will be relatively 
stable since the addition or subtraction of a few clouds will 
not affect the overall shape of that segment of the line. In 
addition, the crossing time at a lower velocity will be longer 
than average. The net result is a stable line core with wings 
that occasionally balloon into shoulders. Shouldering may 
therefore indicate the movement of a few high velocity ex- 
tremely MHD broadened clouds. An observer might be able 
to track individual clouds by subtracting shouldered spectra 
from a mean spectrum or from a pre or post shoulder state 
of a line wing. Based on our modeling of the line profile for 
~ 300 clouds and comparing it with the the amplitude of 
the shouldering observed in NGC 5548 there may be only a 
few hundred clouds in this object (see also Wanders (1997)). 

Current virial mass estimates based on profile width will 
be too big. Energy equipartition of gravity with the magnetic 
wave energy results in an error of a factor of ~ 1.6 2 (~ 2.6) 
in the virial mass if the FWHM of the line is used. This is 
a consequence of Equation 7 (ob ~ v G )- Line profiles are 
wide because of roughly equal contributions from the veloc- 
ity field of the cloud ensemble and the line widths of indi- 
vidual clouds. Additional non-kinematic mechanisms used 
to explain the extreme wings of broad line profiles, e.g. elec- 
tron scattering (Emmering et al. 1992, Bottorff et al. 1997), 
could also be at work though magnetic broadening may ob- 
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viate the need of them to account for the wings of line pro- 
files. The widths of clouds in the extreme wings guarantee 
overlap in one wing with the profiles of clouds in the other 
wing and the line core. The response of the extreme wings 
to short term changes in the continuum (not to be confused 
with the long term kinematic changes discussed above) will 
be smoothed. Thus, for short term variability, due to rapid 
changes in the continuum, our model predicts a variable cen- 
tral part of the profile and a temporally smoothed response 
in the extreme wings. 

An observer's orientation with respect to the magnetic 
field will affect the observed profile if the field is coherent. 
The FWHM of a cloud should be wider in the plane of oscil- 
lation than if viewed along the magnetic field. In addition, 
cloud models need to be modified to include the effects of 
large non-dissipative internal motions on radiation transfer 
within a cloud. 

Finally, dynamical mechanisms for initiating oscilla- 
tions need to be investigated. Spatial wave amplitudes of 

A^\/2n, (12) 

where A is the wavelength, are required to produce trans- 
verse displacement velocities comparable to ob- Since the 
length of a cloud is I « A^/mo x 10 14 cm, where N24 is the 
column density in units of 10 24 cm~ 2 and mo is the cloud 
particle density in units of 10 10 cm~ 3 , the longest wavelength 
supported in a cloud is of the order I « A. Using equation 
(7) gives the period of oscillation, T, is bounded by 

T ^ — ^ 3 x 10 5 — l/^p- s. (13) 

Thus, T is about 4 days or less for fiducial values. We note 
that many continuum engines vary on a time scale of this 
order so radiation pressure could be important. The impor- 
tance of radiation pressure effects in static magnetic out- 
flows is clearly demonstrated in Konigl and Kartje (1994). 
Whether or not oscillations can be initiated by it however is 
still an open question. 



3 CONCLUSIONS 

Magnetic confinement and associated MHD wave broaden- 
ing has broad ramifications. 

• Energy equipartition between non-dissipative waves 
(l/2pa B ) and the magnetic field (B 2 /8ty) seems to be 
an inevitable consequence of magnetically confined gas. 
Line broadening is expected everywhere magnetic fields are 
found. In cases where there is also energy equipartition with 
gravity, the effects of magnetic line broadening will be ex- 
treme. 



• The smallest resolved component in a line may be con- 
siderably broader than its thermal width. If a component is 
part of an ordered flow, then its width can rival that of the 
flow if the magnetic field is in energy equipartition with grav- 
ity. This has important implications for current attempts to 
detect and count individual cloud elements. 

• Line profiles broadened by a magnetic field will be more 
symmetric than expected by predictions of a kinematic 
model that uses thermal broadening only. Thus the param- 
eters describing the shape of a model line profile may need 
to be changed considerably once correct radiative transfer is 
applied. 
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• The spectrum emitted by clouds subject to MHD wave 
broadening needs to be reevaluated since the local velocity 
field will be highly supersonic, not thermal as previously 
assumed. This will have fundamental effects on the spectrum 
since line trapping will be far less severe and continuum 
pumping far more important (Ferland 1999). 

• Virial mass estimates will be too large by a factor ~ 2 if 
the magnetic field is in energy equipartition with gravity. 



